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A recent neutrino experiment at Daya Bay gives superior data of the distribution of the prompt 
energy. In this paper, the energy distribution presented in the experiment is simulated by applying a 
Ganssian-like packet to the neutrino wave function received by the detector. We find that the wave 
packet of neutrinos is expanded during the propagation. As a result, the mixing angle ^13 is more 
difficult to be measured than O 12 and 623 in long baseline experiments. Some other propagation 
properties, such as the time evaluation of the survival probability, the neutrino oscillation and 
the CP violation, are also studied with the employment of the coherent state method. When the 
Gaussian packet width increases, the amplitude of the neutrino oscillation decreases, whereas the 
oscillation period increases gradually. 


I. INTRODUCTION 

Neutrino masses and flavor mixing have attracted 
much attention in the past few years in fundamental 
physics [1, 2] . The mixture is governed by the 3x3 unitary 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix[3- 
5] with three mixing angles (012,013,^23) and a CP- 
violating phase S. Theoretically, the CP violation and 
the neutrino oscillation are experimental observable if 
the mixing angles 0 ij and the squared mass differences 
Am?- = mj — m3 are nonzero [ll. 

Very recently, the nonzero mixing angle 0^ has been 
established for 0i2 and 023 [2, 6, 7] with high accuracies. 
Relatively, 0i3 is more difficult to be determined in exper¬ 
iments than other two angles. G. L. Fogli et al. found 
the strong evidence ( ^ 3 ct level, where a is the stan¬ 
dard deviation) for the nonzero 0i3 in 2011 [8]. Then the 
short-baseline reactor experiments in Daya Bay [9] and 
in RENO [10] improved the confidence level up to 5 ct for 
013 > 0. The Daya Bay experiment found that sin^ 0i3 
0.023±0.003 via the survival probability of the electron 
antineutrino Dg at short distances from the reactor, which 
is given as Pgur = 1 ~ sin^(20i3) sin^(1.267Am^3 • L/E), 
where E is the neutrino energy and L is the travel dis¬ 
tance. The RENO experiment received a very close re¬ 
sult sin^ 013 ~ 0.029±0.006 via a similar setup. However, 
the thoretical investigation for short distance measure¬ 
ment is still incomplete. In this paper, we use the wave 
packet method to study the neutrino oscillation, and find 
that (1) the oscillation is dominated by 0 i 3 in short dis¬ 
tance, (2) the wave packet decay exponentially for large 
distance. Both these results demonstrate that short dis¬ 
tance meansurement is important. 

The accurate determination of the three mixing angles 
is of essential importance to research the CP violation 
in the neutrino sector. Lots of attempts are tried to im¬ 
prove the precision of 0i3. The accuracy of mixing angles 
from oscillation measurement hangs on the availability 
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of a more accurate expression of survival probability, of 
which the implementation needs fewer assumptions, thus 
allowing a flexible experimental arrangement, e.g., the 
source-detector distance is no more a critical prerequi¬ 
site. Such an expression for survival probability should 
be obtained in a strict sense. 

In both the experiments of Daya Bay and RENO, we 
noticed that the surviving probabilities have been given 
with regard to a single energy value of the neutrino, how¬ 
ever the true neutrino beam may have a wide energy dis¬ 
tribution. Eor example, the energy spectrum of the elec¬ 
tron antineutrinos from the reactor at Daya Bay peaks at 
^ 3.0 MeV with a half width of comparable value[9, 10]. 
More details about the mentioned experiments show that 
the neutrino beams utilized have Gaussian-like energy 
distributions (See figure 5 in Ref. [9]). Therefore, a 
more precise and practical method to study the neu¬ 
trino properties presented in the experiments is using 
the Gaussian-like packet to simulate the wave function 
received by the detector. The wave packet approach[ll] 
to the neutrino oscillation has been proposed with differ¬ 
ent motivations[12-15], e.g., to cope with the uncertainty 
principle, but it has not got the deserved attention. In 
this work, we use the Gaussian-like packet to simulate the 
energy distribution of the electron neutrino of Daya Bay 
and RENO, and to study the neutrino oscillation prop¬ 
erties. Eor a given single energy value (say 3.0 Mev) and 
other reliable parameters, we find the electron neutrino 
oscillation period is about A21 = 97.63km, while the fast 
oscillation period is approximately A31 ~ A32 = 3.194km. 
The CP violation in this case has a same oscillation pe¬ 
riod as A 21 . When the energy distribution is taken into 
account, the survival probability of the electron neutrino 
which oscillates along the propagation direction is calcu¬ 
lated, with which we demonstrate that the mixing angle 
013 is more difficult to be measured than 0 i 2 and 623 in 
long-baseline experiments. During the propagation, the 
wave packet of neutrinos gradually disperses with its fluc¬ 
tuation amplitude decreases and the periods simultane¬ 
ously become larger. On the other hand, the similar law 
is found for the oscillation when in terms of the neutrino 
packet width. As initial packet width increases, the os- 
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cillation periods for both the survival probability and the 
CP violation increase, whereas the amplitudes tend to be 
smaller. Numerical results concerning the wave packet 
spreading imply that a sharp neutrino wave packet will 
benefit the detection of the oscillation. 

The present work is organized as follows. In Sec. II, we 
introduce the model Hamiltonian of the one-dimensional 
Dirac equation for massive fermions and the parity opera¬ 
tor, then solve the model Hamiltonian and get the propa¬ 
gation of the Gaussian-like packet wave function. In Sec. 
HI, we calculate some important physical quantities: the 
time evaluation of three-flavor neutrino oscillations, the 
ratio of the density of the electron neutrino, and the os¬ 
cillation of the CP violation of the muon neutrino. A 
summary is given in Sec. IV. 


II. MODEL HAMILTONIAN AND WAVE 
FUNCTION 

In the following context, we assume that the neu¬ 
trino can be described by a wave packet, and the wave 
packet is propagating along the z-diretion. The mo¬ 
mentum components orthogonal to the propagation di¬ 
rection are much smaller than the parallel component 
iiPx), (Py) ^ (Pz)), SO that we can simplify the Dirac 
equation (a ■ p + m/3)'tp = Etjj into a one dimensional 
problem: 


{ctzPz + rnl3)4i = Eijj, (1) 


where the Dirac matrices take the following forms: 


a. 


0 




( 2 ) 



FIG. 1. The electron neutrino oscillation as a function of 
the distance away from the source. The energy of the initial 
electron neutrino is chosen to be the center of the Gaussian 
packet, Aao = S.OMeV. 


Now we investigate the eigenstates of the model Hamil¬ 
tonian. It is easy to check that H commutates with the 
parity operator H = exp (ittN), where N = az + ^ + b^b. 
With some tedious but straightforward derivations, it 
could be proven that the common eigenstates of the 
Hamiltonian and the parity operator can be written as 


|n = +i) = 


|n = -i) = 


^g-f'Vt/2-l-abt _ g-bGt/2-aht^ |Q^ 

( 6 ) 

p_ ^g-bGV2-|-Qh^ _|_ g-hGV2-ab^^ |g^ 

Q f^-b^bC^+ab^ _ g-hGt/Z-afet \ |q^ ’ 


j = x,y,z, m is the mass of neutrino, E is the eigen- 
energy and V' is the wave function with four compo¬ 
nents. One can find that the spin-up and spin-down 
states are decoupled in the one dimensional Dirac equa¬ 
tion (1). Considering the neutrino oscillation problem 
concentrated in this paper is irrelated to the spin degree 
of freedom and without lose of generality, we only take 
into account the spin-up states and then the Dirac equa¬ 
tion can be simplified to be 

f Pz'4’2 + mV'i = E'lpi , , 

\ Pz'4’1 - rml)2 = E^2 ’ ^ 

where ipi , '02 are two components of spin-up states. The 
conjugate operators pz and z can be represented by the 
following bosonic creation and annihilation operators: 


where Tlf, G+, and G_ are arbitrary coefficients, a 
is the eigenvalue of the annihilation operator, and |0) is 
the ground state of the operator h^b with the formula 
b^b\0) = 0. For a given a, energy eigenvalues could be 
found by solving the eigenvalue equation HjH = ± 1 , 0 :) = 
= ±l,a), obtaining 


E^=^^ = ±Vm2VPA^, (8) 

where A = l/-\/2A. The coefficients (unnormalized) cor- 
responding to the positive and the negative eigenstate 
energies are 


F+(_) = l, G+(_) =o/fy2A(m±|£;„|) 


E(a) > 0, 

(9) 


z =^(6-6t) 

(b + b^) 


Pz = 


V2A 


( 4 ) 


where A is a arbitrary number. Consequently, the model 
Hamiltonian (3) could be rewritten as 


E =-^^^{b + b^)+maz. (5) 


F+{-) = -a/ V2A{m + \Ea\) , G+(_) = 1, < 0, 


( 10 ) 

where ±(—) denotes the parity. From Eqs. (9) and (10), 
we can see that the obtained eigenstates are two-fold de¬ 
generate due to the parity symmetry. 

It could be tested that the eigenstates Eqs. (6) and (7) 
with a single vale of o cannot be normalized, so that the 
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FIG. 2. The muon neutrino oscillation as a function of dis¬ 
tance away from the source. The energy of the initial electron 
neutrino is chosen to be the center of the Gaussian packet, 
Aao = S.OMeV. 


FIG. 3. The oscillation of the CP violation of the muon neu¬ 
trino as a function of distance. The energy of the initial elec¬ 
tron neutrino is been chosen to be the center of the Gaussian 
packet, Aao ~ S.OMeV. 


real wave function which is physically significant should 
be a wave packet consisting of eigenstates with differ¬ 
ent values of a. Considering that the neutrino beam 
in the experiments have a Gaussian-like energy distribu¬ 
tion, the neutrino wave functions with definite parities 
are supposed to have the following forms 

|z/(±)) =y’dQ;e-T'(“-“«)'|n = ±l), (II) 

where ap is the center of the Gaussian distribution, 7 
is determined by the half-width, |n = ±1) are given in 
Eqs. ( 6 ) and (7), and F± and G± are chosen to be the 
coefficients of the positive energy solutions in Eq. (9) 
for simplicity. The wave function constructed above is a 
stationary one which is impossible to propagate along 2 ; 
direction. It could be computed that the average prop¬ 
agation velocity of a definite parity Gaussian-like packet 
is exactly zero. For example, the wave function of the 
parity-positive state at time t can be written as 

|zz(+)(t))= = (12) 


Then the expectation values of the momentum operator 
Pz and the position operator z are both zero: 

(z/(+)(t)|p^|i/(+^(t)) = {t)) = 0. (13) 


So that the propagation wave packet must be a mixture 
state of the parity-positive and the parity-negative states. 
Without lose of generality, we choose the following initial 
state whose propagation velocity is maximised: 


= J da 


| 0 ) 

_ a _ 

V2A{Tn+\E^c,)\) 


| 0 ) 


(14) 


In experiments, the neutrinos are found to be of three 
flavours, Vp = (z^e, Vp, Vt)-, when taking part in the weak 


interaction processes. Within the precision of the current 
experimental data, three flavour eigenstates are orthog¬ 
onal with each other. Theoretically, the flavour states 
are linear combinations of the Dirac mass eigenstates 
Vj = {vi,V 2 ,v^)i having masses rrij = (toi, 7712 , 7713 ): 

Vp = (1^) 

where I/pmns is the PMNS lepton flavour mixing matrix 


/I 0 0 \ 

GpMNS = 0 C 23 S23e*‘^=3 

\ 0 C 23 / 


/ Ci3 

0 


f C 12 

si2e*‘^i= 

o\ 

0 

1 

0 • 

-si2e-*‘^i= 

C 12 

0 


0 

Cl3 / 

^ 0 

0 

1/ 


(16) 


in which p = (e,^,r), = cos0y, Sij = sin0y, i,j = 

1,2,3, 9ij are the mixing angles, and (l)ij are the CP 
violation angles. Finally, our initial state of neutrino 
can be rewritten as 


\y,{t = 0)) = N J 


dae ^2 ^PMNsl^i(a))> 

J = 1.2.3 


Wj{a)) = e 2 J jV J 


| 0 ) 

_ a 

y/2A{mj-\-Ea) 



(17) 


where N is the normalization constant. In the ultrarel- 
ativistic limit, the masses of neutrinos are much smaller 
than their energies mj <g; ~ Aa, so we can take 

the following approximation. 


E. 


'j(a) 


= Aa ' 


2 Aa 


V2A{mj + ) 1 + 'm-j/Aa' 


(18) 

(19) 


Accordingly, the normalisation constant and the time 
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suits given in Ref. [9] as Auq = S.OMeV, 7 = 0.056MeV^. 

Firstly, we consider the time evaluation of the proba¬ 
bility of an initial electron neutrino Pe(<a, t) = |' 0 e(<a, 
where the amplitude '0e(<a,t) = {ve\i'{t)). After some te¬ 
dious but straightforward derivations on the basis of the 
Eq. (17), we found 

Pe{.a,t) = iV 2 e- 27 (a-ao)" + ^3 

+ 2 cf^cf 2 Si 2 cos -|- 2 C 12 C 13 S 13 cos (flist) 

+ 2s?2C?3'Si3 cos (^230] > (22) 

and then the total probability corresponding to the 
Gaussian-like distributed electron neutrino beam reads 


FIG. 4. The survival probability of the electron neutrino 
SP{x) = pe{x)/pe{0) as a function of distance for the gen¬ 
eral case (the bine line), the short range approximation 612 = 
023 = 0 (the pink line), and the long range approximation 
^13 = 0 (the yellow line). A Gaussian-like wave function is 
applied to simulate the neutrino beam. 


Pe{t) = J pe{a,t)da = Cl 
+ 2N‘^ j [C 13 C 12 S 12 cos (^ 12 ^) 

+C?2C?3Si3 cos {flist) + Si2C?3S?3 COS (^230] i (23) 


where Ci is a time (distance)-independent constant, 
evaluation of the wave function can be described as, which is given by 


iV2 


o 27 “o/( 47 - 1 ) 


47 — 1 


„2 2 
C 12 C 13 


(1 + mi/AoLnY 


„2 2 

013*12 


+ 


’13 


( 20 ) 


(1 -I- ■ 012 !AaoY (1 + mn/AanYl ’ 

\ve{t)) = ^ j [ci2Ci3e"*-®^(“)Vi(a)) 

+ ci3Si2e*'^i='e“*®2(“)*|i/2(a)) +si3e*'^"^e“*-®=*(“iV3(a))] • 

( 21 ) 


III. PHYSICAL QUANTITIES 

Based on the time evaluation of the wave function 
Eq. (21), we study some interesting physical quanti¬ 
ties. The parameters used in this paper are consistent 
with current measurements: ( 1 ) the mass differences 
Am ?2 = 7.59xl0-5eV2[16], ATOi 2 = 2.32 x 10-W[17], 
Am §3 = Am§ 2 ; ( 2 ) mixing angles sin^ 20 i 2 = 0.861[18], 
sin2 26li3 = 0.092[9], sin^ 26»23 = 1.0[18]; (3) CP phases 
^12 = 4’23 = 0, 013 = 7r/3[19]; (4) Parameters that re¬ 
lated to the Gaussian packet are estimated from the re- 



and the three oscillation frequencies are defined as flij = 
(mf - m]) /(2Aa), {i,j = 1,2,3). 

In practical calculations, the parameter a has the di¬ 
mension eV/A, where rrij <C Aa. We can use Aa as the 
leading order approximation of the energy. In addition, 
the speed of neutrino propagation is chosen to be the 
speed of light. Fig. (1) gives the density evaluation of the 
electron neutrino with the given energy Aao = 3.0MeV 
(the center of the Gaussian-like energy distribution in the 
Daya Bay experiment) as a function of distance. In the 
plot, the slow oscillation period reflects the mismatch of 
the masses difference Amf 2 , while the fast oscillation is 
governed by and ATO 23 . More carefully, the oscil¬ 

lation period is estimated by Ay = dTrAao/Amf^-, which 
gives A 21 = 97.63km and A 31 ~ A 32 = 3.194km. 

Similarly, we can also calculate the oscillation of the 
other two flavor neutrinos. We find that the correspond¬ 
ing probabilities of the muon neutrino and tauon neutrino 
are given by 


Pf,{a, t) =C 2 + 2N'^e [ C?2Cl3*12(* 13*23 “ C 23 ) COs(Qi 2 t) - C? 2 C? 3 Si 3 S ^3 COs(Qi3t) - C?3 Si 2*13*23 COs(Q23l) 

+ C 12 C 23 S 13 S 23 C 33 S 32 COs(Qi 2 t — 012 — 023 + 013) ~ C 23 S 12 S 13 S 23 C 12 S 13 COs(Qi 2 t + 012 + 023 ~ 013) 

~ C 12 C 23 S 12 S 13 S 23 C 13 COs(Qi 3 t — 012 — 023 + 013) + Ci2C23*12*13*23Ci3 COs(Q23l ~ 012 ~ 023 + 013) ] j 

Pr{a, t) =C3 + {ci2C?3*12(*13C23 “ * 23 ) COs(Qi 2 <) - Ci3C23S?3[Ci2 COs(Qi3t) -f S 12 COs(Q23t)] 

+ Ci2Ci3C23Si 2S13*23 COs(Qi2t + 012 + 023 ~ 013) ~ C12C13C23S12S13S23 COs(Qi2t — 012 — 023 + 013) 

+ C12C13C23S12S13S23 COs(r2i3t — 012 — 023 + 013) —Ci2Si 3C23S12S13*23 COs(Q23l ~ 012 — 023 + 013)} , 


(25) 


( 26 ) 
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where C 2 and C 3 are time(distance)-independent constants 


C 2 

C3 


-Ar2p-27(«-ao)^ [0^2 ^2 2 2 , 4 2 2 2 , 2 4 2 2 , 2 2 2 2 

—K L^''12‘'13''23*12 w Ci2Ci3ii3623 ‘'13'’12*13'’23 7^''13'’13*13'’23 

+ 2Ci2Ci3C23Si2Si3S23 C 0 s (^12 + ^(>23 ~ 'Pls) ~ 2ci2C23Si3S23Ci3Si2 C 0 s ((/>12 + (/)23 ~ ^ 13 )] 

=^^e"27(«-“o)" [cf2C?3C^3S?3 + 2c?2C?3S?2S23 + C?3C^3st2S?3 + C?3C^3S?3 
+2ci2C^3C23Si2S13S23(Si2 ~ C 32 ) C0s((/)i2 + 4>23 — ?^'13)] ■ 


(27) 

(28) 


Another important quantity we focus on is the CP 
violation of neutrinos which plays an important role in 
the study of weak interactions. The transition amplitude 
of anti-neutrino, i.e., ^ is obtained by inversing 

the CP violation angle (pij —>■ —ipij- From Eq. (25) we 
can get the ratio of CP violation of the muon neutrino 

Ptj.{a,t) - Pfl{a,t) = 4A^^Cl2C33C23Sl2Sl3S23e“^^^““““^^ 
sin(())i2 -I- (j)23 - 4 >i 3 ) [sin(f2i2t) - sin(f2i3t) -|- sin(f223i)] ■ 

(29) 

The numerical result is shown in Fig. (3). One can find 
that the oscillation period is the same as the period of 
the density oscillations A21 shown in Fig. (1) and (2). 
However the maximum value of CP violation is corre¬ 
sponding to the zero point of the density oscillation, i.e., 
distance ~ 75km. 

The preceding paragraphs are concerned with the 
propagation behavior for neutrinos with the single value 
of the energy. In this case, neutrino quantities fluctuate 
in the space like simple harmonic oscillators with their 
amplitudes and periods keep fixed. Now we concentrate 
on the total effect of neutrinos with the Gaussian-like 
energy distribution. In Fig. (4), we show the survival 
probability of the electron neutrino SP{x) = pe{x)/pe(0) 
as a function of the distance away from the source. The 
Gaussian integration in Eq. (23) is difficult to be carried 
out analytically, therefore it was computed numerically 
by the Romberg method [20]. When the propagation dis¬ 
tance is short (x —?► 0), one can see from Eq. (23) that 
the variation of the survival probability 

dSP(x) 

-;-oc — Ami2 cos^ 013 sin^ 20i2 — Am?3 sin^ 20i3 

xdx 

(30) 

and the term proportional to sin^ 20i3 plays a leading role 
in the oscillation because Amf2 ^ In this case, 

the short range survival probability SPs{x) could be ob¬ 
tained by the approximation of setting 0i2 = 023 = 0. 
On the other hand, when the electron neutrinos oscil¬ 
late to another flavor over a relatively long distance, the 
long range probability SPi{x) could be gained by setting 
013 = 0 due to the fact that sin(0i3) <C sin(0i2). In the 
plot, it could be seen that SPi{x) show good agreement 
with SP{x) when x is large, which indicates the contri¬ 
bution of 013 is small in this situation. This may explain 
why 013 is more difficult to be measured than 0i2 and 
023 in long baseline experiments. On the contrary, it also 



FIG. 5. The survival probability of the electron neutrino 
SP{x) as a function of distance for different widths of the 
Gaussian packet, (a) The long range propagation. The am¬ 
plitudes decrease with the increasing Gaussian packet width, 
while the periods remain unaltered, (b) The short range be¬ 
haviour. The fast oscillations fade out soon for all initial 
packets. 


could be deduced from our results that 0i3 is easier to be 
measured than 0i2 and 023 in short baseline experiments. 
In the Daya Bay reactor neutrino experiment, informa¬ 
tion obtained from a series of detectors, some near the 
reactors and some more than a kilometer farther away, 
are brought together to produce a non-zero value for 0i3 
with a significance of 5.2 standard deviations. Our re¬ 
sults provide a good explanation to their success. 

In Fig. (5) and (6), we study the survival rates of elec- 
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Distance x (km) 


FIG. 6. The oscillations of the CP violation of the muon 
neutrinos as a function of distance for different widths of the 
Gaussian-like packet. 



aia^ 


FIG. 8. The hyper-exponential decay factor g as a, function 
of the widths of the Gaussian-like packet ajaa. 



where /Oe(0), a, and g are fitting parameters. One can 
find that the long range behaver of the oscillation am¬ 
plitude is dominated by the exponential decay factor 
ex.]y[—gx^/[a-\-bx^y\. More importantly, one can find 
in Fig. (8) that the hyper-exponential factor g is al¬ 
most a quadratic fuction of the Gaussian packet width 
ajaQ. These results demonstrate that the signal in short- 
baseline reactor experiments will be more clear. 


Distance x (km) 


FIG. 7. The long range behaver of the electron neutrino os¬ 
cillation for different widths of the Gasuusian packet. 


tron neutrinos and the CP violations of muon neutrinos 
as a function of distance for neutrino beams with different 
widths of Gaussian packets. Contrasting to the prop¬ 
agation of the neutrino with single-value energy whose 
vibration amplitudes are fixed with regard to the dis¬ 
tance X, the Gaussian-like wave packet disperses in the 
propagation with its amplitude diminishes and its period 
increases. In the meantime, when the source beam has a 
larger packet width, the fluctuation amplitudes also low¬ 
ered its value whereas the period increases. The results 
indicate that a wide energy distribution of the initial elec¬ 
tron neutrino obstructs the detection of neutrino oscilla¬ 
tions. In order to investigate the dependance of neutrino 
oscillation on the distance x, we plot the amplitude of 
Pe{x) versus x for different Gaussian packet widths in 
Fig. (7). In order to obtain the quantitative relationship 
between the oscillation amplitudes and the widths of the 
Gaussian packet, we fit Fig. (7) with the decay formula 
of Gaussian packet: 


Pe{x) 


PejO) 
\/l + ax^ 


(31) 


In this paper, we investigate the neutrino propagation 
properties by using the Gaussian-like packet to simulate 
the energy distribution of the neutrino wave function re¬ 
ceived by the detector. The time evaluation of the sur¬ 
vival probability, the neutrino oscillation, and the CP vi¬ 
olation are calculated with the application of the coherent 
state method. Our results indicate that the mixing angle 
013 is more easier to be measured in short baseline ex¬ 
periments than 012 and 023 , with which we evidenced the 
success of the Daya Bay experiment in the determination 
of 013. During the propagation, the wave packet diffuses 
to a wider range and the neutrino oscillation fades. Fur¬ 
thermore, the oscillation amplitude decreases with the 
increasing Gaussian-like packet width, whereas the os¬ 
cillation period becomes larger. This suggests that it is 
more easier to detect neutrino oscillations when utilise a 
sharper distribution of neutrino energy. 
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